Unitary evolution and uniqueness of the Fock quantization in flat cosmologies with 

compact spatial sections 



(N 
(N 



Laura Castello Gomar 

Facultad de Ciencias Fisicas, Universidad Complutense de Madrid, Ciudad Universitaria, 28040 Madrid, Spaii^ 

Jeronimo Cortez 

Departamento de Fisica, Facultad de Ciencias, Universidad Nacional Autonoma de Mexico, Mexico D.F. 04510, Mexic^ 

Daniel Martin-de Bla^ and Guillermo A. Mena Marugan 
Instituto de Estructura de la Materia, CSIC, Serrano 121, 28006 Madrid, Spai'^ 

- - Jose M. Velhinho 

. Departamento de Fisica, Faculdade de Ciencias, Universidade da Beira Interior, 

^ ■ R. Marques D'Avila e Bolama, 6201-001 Covilha, Portugal^ 

, We study the Fock quantization of scalar fields with a time dependent mass in cosmological sce- 

O ' narios with flat compact spatial sections. This framework describes physically interesting situations 

' like, e.g., cosmological perturbations in flat Friedmann- Robertson- Walker spacetimes, generally in- 

I eluding a suitable scaling of them by a background function. We prove that the requirements 

of vacuum invariance under the spatial isometrics and of a unitary quantum dynamics select (a) a 
unique canonical pair of field variables among all those related by time dependent canonical transfor- 
mations which scale the field configurations, and (b) a unique Fock representation for the canonical 
commutation relations of this pair of variables. Though the proof is generalizable to other compact 
' spatial topologies in three or less dimensions, we focus on the case of the three-torus owing to its 

O^' relevance in cosmology, paying a especial attention to the role played by the spatial isometries in 

^ , the determination of the representation. 

1 1 

> . I. INTRODUCTION 

l> ■ 
■ 

\^ The construction of a quantum theory to describe a classical system is a process plagued with many ambiguities. 

I Generically, the correspondence is not one-to-one. What is more important, the quantum physics depends on the 
' choices made in the steps affected by those ambiguities. The specification of a unique quantization must then either 
I ^ rest on the confrontation of the predictions with experiments, or be achieved by appealing to additional criteria, usually 
related to symmetries of the system or to the behavior of the quantum states. The problem is especially relevant 
7-H in cosmology, both because the windows for quantum effects in cosmological observations are certainly narrow (if 
. any), and because one cannot really select the best candidate for a quantum model of the universe by performing an 
. ^ ' indefinite number of repeated measurements (in copies of the same state) , since we can only observe the universe in 
I which we live. 

5h • Even if one identifies a classical system with a specific algebra of observables, attained by first selecting a set 
. y. of appropriate variables for the system and constructing out of them the algebra of functions which we declare of 
interest, its quantization is still not unique. In general, there will exist representations of that algebra which are not 
equivalent (that is, which are not related by a unitary transformation). This problem appears not just for generic 
systems with complicated phase spaces or interactions, but also for the simplest systems with a linear behavior. In 
the case of standard quantum mechanics, when the number of degrees of freedom is finite, the ambiguities in the 
representation can be removed by imposing additional requirements. For the Weyl algebra corresponding to the 
free particle, one usually demands that the representation be irreducible, unitary and strongly continuous. Then, 
the Stone- von Neumann theorem guarantees that the representation is unique (up to unitary transformations) [5]. 
A similar result does not exist when one considers systems with fields, which have an infinite number of degrees of 
freedom. The ambiguities persist even if one takes advantage of the linear structures and restricts the consideration to 
representations of the Fock type, where a concept of particle and vacuum are available (either with a genuine physical 



* llaucaste@estumail.ucm.esj 

t jacq@ciencias.unam.mx 

^ daniel.martin@iem.cfma c.csic.esl 

5 mena@iem.cfmac.csic.es^ 

^ |jvelhi@ubi.pt| 



2 



interpretation or just as mathematical entities). The possible choices of Fock representation which are physically 
different are still infinite, and correspond to non-equivalent choices of a vacuum 

For fields propagating in highly symmetric spacetimes, the symmetries of the background can be employed as criteria 
to determine the quantization [2tJ|, asking e.g. that all the quantum structures incorporate those symmetries. For 
instance, this happens in Minkowski spacetime, where the vacuum can be selected by demanding Poincare invariance 
0. But, for more general spacetimes, no generic criterion exists that specifies the Fock representation. In generic 
situations, there is not enough spacetime symmetry to fix it. Therefore, supplementary or alternate requirements are 
needed in order to pick out a unique equivalence class of Fock representations. Frequently, the situation found in 
cosmology is that the spatial sections still present a high degree of symmetry (at least in a certain approximation) , 
but the stationarity is lost owing to the universe expansion (or contraction). In this framework, it seems natural 
to adhere to the remaining spatial symmetries, demanding that they are naturally included in the quantum theory, 
and replace the criterion of time symmetry with the closest possible one, namely, with the requirement of a unitary 
evolution. Unitarity will guarantee a standard probabilistic interpretation in the quantum theory, though the loss 
of time symmetry will make the vacuum change dynamically. Actually, the combined criterion of invariance under 
the spatial symmetries and a unitary evolution has been successfully employed recently in the selection of a unique 
quantization for cosmological models. 

Apart from the choice of representation, there is an additional ambiguity that arises naturally in the construction 
of a Fock description for fields in cosmological scenarios. The non-stationarity of the spacetime leads to the obvious 
possibility of absorbing part of the time dependence of the field via its scaling by a function that depends only 
on the background. This is so irrespective of whether the spacetime in which the propagation takes place is a 
true physical background [sj, Q , an effective background (emerging from an effective description of the system, for 
instance incorporating some quantum modifications at an effective level 0421]) or an auxiliary background (like, e.g., 
in dimensional reductions of gravitational systems in General Relativity, using the presence of Killing symmetries 
p^|-[l2!|). Scalings of this type are found in many circumstances when considering fields or cosmological perturbations 
around Friedmann-Robertson- Walker (FRW) spacetimes, for instance, like in the case of test fields, or scalar and 
tensor perturbations, including their description in terms of gauge invariants, as with the Mukhanov-Sasaki variables 
[isj . These scalings render the field equations in the form of those corresponding to a Klein- Gordon (KG) field in a 
static (auxiliary) spacetime. The time dependence still shows up via the appearance of a varying mass term. Besides, 
usually there is no friction term in the scaled equation, or one can neglect it when considering short scales in the 
field behavior. The scaling of the field configuration can always be regarded as part of a time dependent canonical 
transformation in the system, which is natural to consider as a local and linear transformation in order not to spoil 
this type of properties in the field system. In such a transformation, the field momentum gets the inverse scaling 
of that of the field configuration [to maintain the canonical commutation relations (CCR's)] and, in addition, may 
admit an extra contribution linear in the field configuration, multiplied by a function of time. Any of these canonical 
transformations leads to a new canonical pair for the field, but also alters the dynamics, since the change is time 
dependent. The quantum description can be made by adopting any of these canonical pairs for the field, introducing 
a new kind of ambiguity with infinitely many possibilities. 

As we have briefly commented, the criterion of invariance under the spatial symmetries of the classical spacetime, 
together with the unitary implementation of the dynamics in the quantum realm, has been employed in cosmological 
systems to remove the ambiguity in the Fock quantization. Indeed, this criterion has shown useful not only to select a 
Fock representation, but also to determine a unique canonical pair for the field among all those that can be reached by 
means of linear canonical transformations which depend on time. Note that, since these transformations modify the 
dynamics, the fact that the evolution is unitary is intimately related to the canonical pair adopted. The criterion was 
first applied to the case of the Gowdy models [Tl|, [l^, [l3 - [T8| : reductions of General Relativity with two spatial Killing 
vectors and compact spatial sections. They describe gravitational waves propagating in cosmological spacetimes with 
compact universes. In the case of linearly polarized waves, the system admits a description in terms of a KG field 
(with time dependent mass) in an auxiliary, dimensionally reduced stationary spacetime. The spatial sections of 
this auxiliary background can be either isomorphic to the circle, S^, or to the sphere, S"^, depending on the spatial 
topology of the Gowdy model (which can be that of a three-torus, of a three-sphere, or of a three-handle, x S^). 
The conclusions about the validity of the proposed criterion for the choice of a unique Fock quantization in the Gowdy 
models have been lately extended to the case of backgrounds with spatial sections isomorphic to d-dimensional spheres, 
with d < 3 [T9l - [2^ , and even more recently to a general compact topology in three or less spatial dimensions [2^ [IJ] . 
Although this generalization proves that the uniqueness holds in any case with compact spatial sections, we will focus 
here on the case when these sections have the topology of a three-torus. The relevance of this case is clear since it 
describes flat universes, which is precisely the favored scenario for the universe according to observations. Moreover, 
the discussion for fields in spacetimes with generic compact topology [H, HJ] is obscured by the fact that there may be 
no clear geometric interpretation of the considered symmetry group, in the most general case. This interpretation is 
neat in the case of the three-torus. Besides, there exist certain peculiarities which motivate the mathematical interest 
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in the analysis of the three-torus. Namely, since the group of isometries of the three-torus is an Abelian compact 
group (as in the case its irreducible (unitary) representations are one-dimensional and defined over complex 

vector spaces [l^. This introduces certain subtleties in the characterization of the Fock representations which are 
compatible with those symmetries. These subtleties arise because the complex representations of the symmetry group 
must be combined suitably so that, at the end of the day, one deals exclusively with real scalar fields. In order to fully 
take into account this issue, and show explicitly how those complex representations combine, a departure from the 
general approach followed for other topologies in the literature is adopted [l^ [2l| . Explaining the characterization of 
the symmetric representations will be one of the main goals of this work. 

The paper is organized as follows. In Sec. |TT]we start by considering a KG field with time dependent mass in a static 
spacetime with flat spatial sections of three-torus topology. We briefly present the classical system and the standard 
procedure to introduce a Fock quantization of its associated phase space. Then, in Sec. Illll we prove that our combined 
criterion of a) invariance under the isometries of the three-torus and b) unitary implementation of the dynamics, selects 
a unique Fock representation, up to unitary transformations. This unitary class of representations includes the one 
which would be naturally adopted if the field had vanishing mass. Our proof contains a detailed discussion of the 
characterization of the representations which are invariant under the symmetries of the three-torus, with a careful 
treatment of the symmetry transformations and the consistency conditions coming from the reality of the field. We 
go beyond this result of uniqueness in Sec. IIVI where we analyze time dependent canonical transformations arising 
from a scaling of the field configuration, and demonstrate that only one of all those transformations is compatible 
with our requirements of spatial symmetry invariance and unitary evolution: the trivial one. Finally, we present our 
conclusions in Sec. |Vl The discussion in this work follows a line of arguments similar to that presented in Ref. f2(| . 
to which we refer the reader for further details about the analysis. 

II. KLEIN-GORDON FIELD WITH TIME DEPENDENT MASS 
A. The classical model 

Before proving that the criterion that we put forward indeed selects a unique class of unitarily equivalent Fock 
representations for a linear scalar field with time dependent mass in a fiat spacetime with compact sections, let us 
start by describing the classical set-up of our theory. 

We consider a real scalar field ip defined on a flat spacetime whose spatial sections have the topology of a three-torus, 
T^. These sections are equipped with the standard spatial metric of the three-torus, hij{i,j = 1,2,3). The field is 
subject to a linear equation of KG type: 

(f ~ Aif + s{t)ip ^ 0, (1) 

where the dot stands for the time derivative, A is the Laplace-Beltrami (LB) operator associated with the three-torus 
metric /lij, and s{t) can be interpreted as a time dependent mass. In principle, s{t) can be any time function, and 
only later in the discussion we will impose on it some extremely mild conditions about its derivatives. On the other 
hand, the time domain in which the KG field is defined can be any arbitrary connected real interval I C M. We do 
not impose that I be the real line, not even that it be unbounded. No specific form for the interval I is assumed. 
This is important for the applications of our results to situations in which the field description is effective, since the 
validity of the effective spacetime geometry can be restricted just to a certain finite, bounded time interval. More 
generally, the domain I might even be just a union of connected components. In such cases, the restriction to one of 
the components is sufficient to achieve the same uniqueness results. 

This equation can be obtained in very different cosmological models. An important class of systems for which it has 
a major relevance are scalar fields promgating in non-stationary, cosmological spacetimes, like e.g. the case of matter 
fields in infiationary backgrounds [a, [23, [23] ■ Another class are cosmological perturbations around FRW spacetimes. 
According to our comments in Sec. [II a suitable time dependent scaling leads the field equations of those systems to 
a KG equation of the above type [2ir[2^. In this cosmological context, the case of flat spatial sections that we study 
here is the most interesting one owing to its potential applications to describe situations in the observed universe. 

Clearly, this dynamical equation is invariant under the group of isometries of the three-torus, since the LB operator 
is defined in terms of the standard metric on T'^. To incorporate and analyze the role of these isometries, we consider 
the composition of rigid rotations in each of the periodic spatial directions 9i that diagonalize the metric hij , 

: 0, -> 0, + a„ \fa,eS\ (2) 

Here, ai is the angle parameter that provides the rotation in the direction i {i = 1,2,3). We will call T^, with 
d = (ai,a2,Q;3), the transformation obtained by composing the corresponding rotations. 
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The canonical phase space T of the field system is obtained from the Cauchy data at a reference time io G ^, namely 
{{ip,Pip)} = {(<P|t(,,</5|t(,)}, equipped with a symplectic structure ft that amounts to the canonical Poisson brackets 

{(p{6),P^{6')} = S^{6 - 9'), with S^{6) being the Dirac delta on and 9 the spatial point with coordinates 9i. Note 
that we assume that the Hamilton equation for the field momentum is = tp (29j . 

As a result of the periodicity on the spatial coordinates 0^, one can decompose the field ip (and its momentum) in an 
expansion in Fourier modes. Note that these modes are eigenfunctions of the LB operator. The field decomposition 
using these complex eigenfunctions takes the form 

^{t. 0) - J2 1™ W exp{«(m • 9)}, (3) 

where m is the tuple of integers (mi,m2,rn3) (i.e., e Z, for i — 1,2,3), and we have introduced the notation 
rh-9 — m.i9i. Recalling that the field is real, we find that the coefficients of the expansion are subject to the reality 
conditions 



q-m(i) = [qm(i)]*, (4) 

where the symbol * denotes complex conjugation. 

This complex decomposition is well adapted to the three-torus symmetries. The Fourier modes are eigenfunctions 
of all the transformations Tg. Moreover, each tuple m provides a different, inequivalent irreducible representation 
of the isometry group. In particular, we straightforwardly see that such irreducible complex representations are 
one-dimensional, as it corresponds to the case under consideration, with symmetries that form an Abelian compact 
group. The disadvantage of using this complex Fourier decomposition, nonetheless, is that we have to deal with the 
complications posed by the reality conditions. In order to avoid these complications, and clarify how the irreducible 
representations over complex vector spaces combine in the case of real fields, it is then advisable to adopt an alternate 
decomposition in terms of real eigenfunctions of the LB operator, namely, in terms of Fourier modes corresponding 
to sines and cosines: 

^{t, 9) = T^^Y^lquit) cos(n • 9) + Xfi{t) sin(n • 9)]. (5) 

n 

Here, we have defined the tuple n = (711,712,71.3), with Ui £ 1^ {i — 1,2,3). Note that we have changed the notation 
for the labels of the Fourier modes and coefficients from 777, to 77, to facilitate the distinction between the complex 
and the real formulations. The sum in the above expansion contains only tuples 77 with components ti^ that are all 
non-negative or at most contain one negative integer. All different tuples satisfying this restriction are to be summed 
over (and just once each). We indicate the restriction in the sum with a tilde, instead of making it explicit, something 
which would complicate the notation in excess. Furthermore, in the rest of our discussion we will ignore the zero 
mode, 77 = (0,0,0). As we can anticipate, the unitary evolution and the uniqueness of the representation do not 
depend on the removal of a finite number of degrees of freedom. This mode can always be quantized separately, 
including the possibility of employing non- standard methods in its quantum mechanical description (like, e.g., using 
loop quantization methods) . Its exclusion does not alter the field properties of the system, on which we concentrate 
our attention. 

Since we have expanded the field in eigenfunctions of the LB operator, the different modes are dynamically decou- 
pled. For the cosine modes, we obtain 

qri + [ujI + s{t)]qri = 0, (6) 

and similarly for the sine modes, Xft. Thus, the dynamical equations only depend on the corresponding LB eigenvalue 
of the mode, — w^, which is given by = ^ ■ nf. The subindex 77, introduced as a label in w„, is taken as a positive 
integer which designates the order of these eigenvalues. That is, w„ < aj„' if 71 < 77'. Note that, given the compactness 
of the three-torus, the LB eigenvalues form indeed a sequence {a;„}, which is certainly unbounded. It is also worth 
noticing that the degeneracy (/„ of each eigenspace of the LB operator presents a complicated dependence on the label 
77, because of accidental degeneracy: one can find different tuples 77 (others than those related by a flip of sign in one 
of the components, or by permutations of the components) which lead to the same eigenvalue. For the sake of an 
example: the tuples (2, 2, 1) and (3, 0, 0) correspond to the same eigenvalue — 9. 

Finally, we decompose the field momentum in the same way as we have explained for the field configuration. 
Then, the momentum coefficients for the sine and cosine contributions in the real modes expansion, which we will call 
Pf{ and T/fj, respectively, satisfy the dynamical equations pn = QH: Si,iid = i^. The non- vanishing Poisson brackets 
in terms of these coefficients are {qfi,Pn'} = {xft^yfi'} = Snfi'- 
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B. Fock quantization 



The Hilbert space of the system over which we construct the quantum theory is the direct sum of the symmetric 
tensor products of the one-particle Hilbert space Ho- The Fock space, then, is determined by the one-particle space. 
For the construction of the latter, the only ingredient needed is a complex structure J 0]. We recall that a complex 
structure is a real map on phase space whose square is minus the identity and that preserves the symplectic structure 
f2 In addition, we demand the complex and the symplectic structures to be compatible in the sense that their 
composition, r2( J-, ■) must be a positive bilinear form. The sector of positive frequency (complex) solutions is obtained 
with the projector ^(1 — iJ), where 1 is the identity. This sector is then completed into the one-particle Hilbert space 
using the norm provided by ^[rt{J-, •) — •)]• Therefore, the complex structure contains all the information that 
is physically relevant to determine the different possible Fock representations. This can be rephrased by saying that 
a complex structure, together with the symplectic structure, defines a state, which is usually called the vacuum, and 
then a representation of the OCR's. 

In order to discuss the Fock quantization of the scalar field, we introduce a particular, well-known complex structure 
which we will employ as a starting point in our analysis. The chosen complex structure, Jq, is the one which would 
be naturally related to a free, massless KG field. Hence, it is determined entirely by the three-torus metric, via 
the corresponding LB operator. In particular, this property guarantees that the complex structure Jq, and thus the 
vacuum selected by it, is invariant under the isometry group of the three-torus. In terms of the LB operator, we can 
define Jq using a decomposition in eigenmodes, introducing the annihilationlike variables 

an = J — {uj,iqn + ipn), an = ^ — {uJnXn + iyn)- (7) 

The corresponding creationlike variables are given by the complex conjugates at and 5,1. These new variables provide 
a complete set of coordinates on phase space. The action of Jq on this new basis is defined to be diagonal, in the 
standard form 



Jo{an) = ian, 
Man) = -iah 



Jo{an) = ian, 
Jo{an) = -ia'fi. 



(8) 
(9) 



The evolution of these variables from a fixed initial time to any another time i G I can be expressed as a linear 
transformation U (since the field equations are linear). Given that different modes decouple in the dynamics, the 
evolution is actually block diagonal, with 2x2 blocks C/,-j, one for each pair of annihilation and creationlike variables. 
Furthermore, since the dynamical equations of the modes depend only on the eigenspace of the LB operator under 
consideration (labeled by n), the same happens with the mentioned blocks, which we will therefore designate by J7„. 
Therefore, it is straightforward to conclude that the evolution can be described in the form 

( aft{t)\ _ ( afi{to)\ Tj _ ( o:„{t,to) l3n{t,to)\ , , 

I <it) y " I ^nito) j ' ^" - I3:it^ to) to) ) ' ^'"^ 

and similarly for (a,-j,al). Finally, since the dynamical evolution is a symplectic transformation, the alpha and beta 
functions appearing in this matrix expression must satisfy the relation 

\an{tM)? -\f3n{tM)? ^l, Vie I, (11) 



for all values of the eigenvalue label n. 

In the next section we will show that the Fock representation determined by Jq admits a unitary implementation 
of the evolution and, furthermore, that this is the only representation with that property (up to unitary equiva- 
lence) among all those that are invariant under the isometry group of the three-torus. We recall that, in general, 
a linear canonical transformation U can be implemented quantum mechanically as a unitary transformation in the 
representation determined by a complex structure J if and only if the antilinear part of the transformation, given 
by Uj = ^{U -\- JUJ), is a Hilbert-Schmidt operator fs^], namely, that the trace of UjUj is finite (here, the dagger 
denotes the adjoint operator). It is possible to rephrase this condition as the requirement that the antilinear coeffi- 
cients of the considered transformation (usually called the beta Bogoliubov coefficients) be square summable (that is, 
that their squared norms have a convergent sum) . Another equivalent way of stating this condition is to demand that 
the image of the vacuum state under the transformation U possess a finite number of "particles" , using the particle 
concept associated to the original vacuum. We will use this condition for unitary implementability (in any of its 
versions) in the rest of our discussion. 
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III. UNIQUENESS OF THE QUANTIZATION 



We will now show that the Fock representation determined by the complex structure Jq leads to a unitary quantum 
evolution, even if the field has in fact a time dependent mass. We will also characterize the most general complex 
structure that is invariant under the symmetry group formed by the transformations T5. We will see that they are all 
related by a specific family of symplectic transformations. Using that characterization, we will prove the uniqueness of 
the invariant complex structure (up to unitary equivalence) under the requirement that the dynamics admit a unitary 
implementation. To avoid repeating parts of the demonstration which follow the line of arguments presented in the 
literature for the three-sphere [l9| . we will review the main steps of the proof and concentrate our attention just on 
the aspects that are specific of the three-torus. 



According to our comments above, the evolution U is implementable as a unitary transformation on the Fock space 
determined by Jq if and only if its beta coefficients /3„(t, io) [appearing in Eq. (jlOp ] are square summable. Taking 
into account the degeneracy gn of the eigenspaces of the LB operator, the necessary and sufficient condition is the 
finiteness of the sum gn|/3n(i, io)P for all possible values of time tel. Physically, this amounts to a finite particle 
production during evolution. Clearly, the summability of the sequence of beta functions (in squared norm) depends on 
the asymptotic behavior of /?„ (and gn) in the ultraviolet, namely, in the limit of infinitely large eigenvalues ujn- The 
asymptotic analysis of the Bogoliubov coefficients for a KG field with time dependent mass in a stationary spacetime 
was carried out in Ref. [19j], and the result applies in particular to the case discussed here, when the compact spatial 
sections are isomorphic to a three-torus. The analysis leads to the conclusion that 



The symbol O indicates the asymptotic order. The only hypothesis about the field that is employed to deduce this 
behavior (but that it is not even necessary for its validity) is that the mass function s{t) possesses a first derivative 
which is integrable in every compact subinterval of I. Using the asymptotics, it is straightforward to see that the 
condition for a unitary implementation of the evolution is equivalent to the summability of the sequence formed by 
gn/^jJn- To check whether this summability holds, one has to study how the degeneracy g„ changes with n in the limit 
when this label gets infinitely large. This variation of gn is quite involved owing to the accidental degeneracy that we 
have already pointed out. The exact dependence of the degeneracy with n cannot be given explicitly. Nonetheless, 
for our discussion, we only need to compute the asymptotic behavior of gn, a task which can be done in a relatively 
simple way as follows. 

The values can be understood as the norm of the vector ft provided by the tuple that labels the modes. In 
principle, n is restricted to the octants in with none or at most one negative component. However, there exist two 
modes for each value of n: the sine and cosine modes. We can assign the two modes to the couple of vectors (n, — n). 
It is then clear that we can make correspond modes to vectors with integer components (except the zero-mode, that 
has been excluded). Let us call the number of modes whose eigenvalue uJn is in the interval {N,N + 1], with N 
a natural number. Since is strictly decreasing with n, we have that the sum J2ni9n/'^t)j whose convergence 

we want to check, is always equal or smaller than J^wi-^N /N'^). Geometrically, the value of Dn is the number of 
vertices of the cubic lattice with step equal to one that are contained between the sphere of radius + 1 (including its 
surface) and the sphere of radius N. Therefore, Dn increases with N like A^. It is then straightforward to see that 
the sum of D^/N^ is finite, and a fortiori that of gn/uJn- Thus, the Fock representation determined by Jq, naturally 
associated with the case of a massless field, provides in fact a unitary implementation of the dynamics even in the 
presence of a time dependent mass. 



We will now characterize the most general complex structure that is invariant under the group of symmetries 
formed by the transformations Ts, corresponding to the isometrics of the three-torus. For simplicity, we will just 
call invariant such complex structures. To reach this characterization, we will follow a procedure which differs from 
that presented in the literature for other compact spatial topologies, adopting our analysis to the peculiarities of the 
three-torus. As we have remarked, in our case, the isometry group is Abelian, and hence its irreducible representations 
are all one-dimensional over complex vector spaces. On the other hand, since the scalar field studied is real, these 



A. Unitary evolution in the massless representation 




(12) 



B. Characterization of the invariant complex structures 
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representations must be combined suitably. To clarify this interplay, we will rather consider real representations from 
the start, adapted to the decomposition of the field in sine and cosine modes. 

It is also worth noticing that, owing to the commented accidental degeneracy, the eigenspaces of the LB operator 
do not provide irreducible representations of the isometry group, even if the operator commutes with the isometrics 
because it is constructed out of the three-torus metric. Again, this situation is novel compared to that found in the 
literature, e.g. for d-spheres (T9l - [2^ . Hence, the case of the three-torus calls for a detailed analysis, that we present 
in the rest of this subsection. 

The action of the three-torus symmetries T5 on the real sine and cosine modes [see Eq. ([5])] is easily derived from 
its active transformation of the field. One gets 



q'fi \ _ f cos{fi ■ a) ~ sin(n • a) \ / qrx 
sin(n • a) cos(n - a) I \ Xfi 



(13) 



Similar equations are obtained for the modes of the field momentum (p^j, y^). We notice that these transformations 
only mix modes in pairs, just the sine and cosine modes with the same label n. Furthermore, the action on each of 
these pairs is different. The sine and cosine modes that get mixed belong to the same LB eigenspace and hence have 
the same dynamics. Besides, as we anticipated, not all modes in the same eigenspace get mixed under the action of 
the symmetry group, owing to the accidental degeneracy. This fact complicates the characterization of the invariant 
complex structures and leads to a situation which is similar to that encountered for the topology US] : though 
in higher dimensions. Examining the action of the transformations Ts given above, it is not difficult to realize that 
the most general invariant complex structure must be block diagonal in the label rt, namely, 

J = ^'Jn, (14) 

n 

where each complex structure corresponds to a 4 x 4 block, associated with the sine and cosine modes determined 
by n for the field configuration and momentum. In the direct sum, we have used the same kind of notation introduced 
in Eq. ©. 

Let us express the blocks Jn in terms of smaller 2x2 blocks in the phase (sub)space basis formed by {qn, xn^Pn, TJa)'- 

From the condition that J be an invariant complex structure, so that TZ^JT^ — J for all transformations Tg, one 
concludes after a straightforward computation that every 2x2 block {Qa} — {Aji^ Bfi,Cfi, Dn} must commute with 
all the rotation matrices of the form 

(16) 



sin(n • a) cos(n • a) 
that is 

R^\d)QnRn{3)^Qn- (17) 

It then follows that every 2x2 block must have a diagonal part proportional to the identity and a skew-symmetric 
non-diagonal part, namely, 

where Q^^^ and Q^^' are arbitrary real numbers. 

We still have to impose the condition that the invariant complex structure be compatible with the symplectic 
structure, so that their combination ri( J-, •) provides a positive definite bilinear map on phase space. In terms of the 
blocks Jfj, this condition implies that Jjils must be a positive definite symmetric matrix, where is the transpose 
of Jfi and the blocks of the symplectic structure are 

\ 12x2 Ll2x2 / 

Here, 02x2 is the zero matrix in two dimensions and l2x2 is the identity matrix. In order to satisfy this requirement, 
the 2x2 blocks of Ja given by Bfi and Ca must be symmetric matrices of negative and positive definite type. 
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respectively, whereas the two other blocks must satisfy that Afi = — D^. Together with the condition of invariance 
discussed above, we then conclude that the blocks J^j must have the form 

where Bf^ and must be proportional to the identity, with a non-positive and a non-negative constant of propor- 
tionality, respectively. That is, 

^(2) ^ ^,(2) ^ Q ^{1) ^ Q ^{1} ^ Q .21) 

In addition, we must also impose the remaining conditions that the square of the complex structure be minus the 
identity, = —1, and that it leave invariant the symplectic structure, namely f2(J-, J-) = ■). These requirements 
lead in turn to the following conditions on the blocks J^: 

From the former of these restrictions we obtain the following equations for the matrix elements: 

-{^Y^l^Y-B^C^^^^ 4^^f =0. (23) 
On the other hand, the second restriction leads to the equations 

^B^=^^ 4'4^^=0. (25) 

Summing the first equality of Eq. ((23)) and Eq. (p4|) . we get that the matrix element mush vanish. Therefore, 
every block J,-j of an invariant complex structure (compatible with the symplectic form) must have the form 



.^_( 4^^12X2 4^'l2x2 \ /r,fi^ 



where 



B^'<^, CW>0, A^=±^-l-B^C^. (27) 



Note that B^}^C^}^ < ^ 1, since A^}^ must be a real number. 

Once we have deduced the general expression of the invariant complex structures in the basis {(qa, Xfi,pfi, ya)}, it is 
straightforward to write them in the basis {{af{,a^,df{,dV)} of annihilation and creationlike variables defined by the 
complex structure Jq, basis in which we can easily compare them. The considered change of basis can be obtained 
from Eqs. ([7]). In this way, one can check that the blocks of any admissible invariant complex structure must be 
block diagonal, with 2x2 blocks that coincide by pairs. More explicitly, 

Jn = ( ) ^ (28) 



02x2 >Jn 

with 

^"~M r(1) , r-(l), ,2 , ,-4(1), , r(1) ,2 ■ ^^^1 

Actually, every invariant complex structure J (that is compatible with the symplectic structure) can be obtained 
from Jo by means of a symplectic transformation, K, namely J — KJoK~^ uM- Since, in the considered basis 
of annihilation and creationlike variables, Jq is diagonal with blocks of the form (Jo)n — diag{i, — i, i, — i}, the 
transformation K can also be taken block diagonal, with 4x4 blocks of the type 
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Here, Kfi and are complex numbers which play the role of alpha and beta Bogoliubov coefficients for the transfor- 
mation K. In particular, they satisfy the symplectomorphism condition 

\i^n\^ - |A,ip = 1 Vn. (31) 

Finally, the relation of the coefficients Kfi and A^; with the matrix elements ^^j^"*, b'^\ and cl^'* is given by 

2|^„-|' = l-4'^+4'^^^, (32) 

Notice that the phase of Kft can be chosen freely. For instance one can choose it so that this Bogoliubov coefficient 
be non-negative. Actually, one can see that this choice does not affect the rest of our considerations. 



C. Uniqueness of the invariant representation with unitary dynamics 

To conclude our proof, showing the validity of our criterion to pick out a unique Fock representation, we still have 
to consider the unitary implementation of the dynamics and demonstrate that this restricts the admissible invariant 
complex structures to only one class of unitary equivalence. We will use again the fact that all the possible invariant 
complex structures J are related with the complex structure Jo by means of a symplectomorphism K of the form given 
above, which can be understood just as a change of annihilation and creationlike variables. It follows that the unitary 
implementation of the evolution C/ in the representation determined by J amounts to the unitary implementation 
of K~^UK with respect to Jq [11]. The beta coefficients of K^^UK can be viewed as the antilinear coefficients of 
the Bogoliubov transformation determined by the dynamics, expressed in terms of the annihilation and creationlike 
operators selected by J, instead of by Jq. A trivial computation shows that these new beta coefficients, that we will 
call pi, take the following expression in terms of the original Bogoliubov coefficients: 

PiitM) = {KnfMt,to) - (Xnff3*{t,to) + 2iKt\n3[an{t,to)]. (34) 

Here, denotes the imaginary part. It is worth remarking that these beta coefficients (or rather beta functions, 
taking into account the time variation) depend now not just on n, the label of the eigenspaces of the LB operator, 
but rather on n, which is the label of the sine and cosine modes and, as a consequence, of the "irreducible" reai 
representations of the symmetry group of the three-torus. Therefore, if, according to our criterion, we restrict our 
discussion exclusively to invariant complex structures J which allow a unitary quantum evolution, the above beta 
functions will have to be square summable (over n) at all values of time t £l. So, we will assume that this is the case 
from now on. 

The analysis of this summability can be made along a line of arguments similar to that presented in Ref. (l9j . Using 
that the beta functions Pnit, to) corresponding to Jq are square summable and that jK^j > 1, one can prove that the 
summability (in square norm) of I3i{t, to) implies the same property for the set formed by 3[a„(i, to)]AH/K^. Recalling 
then the asymptotic behavior of an{t,to) and calling z,-j = Ah/kI, we arrive at the conclusion that the quantities 



Z0 sm 



w„(i - to) + 



(35) 



form a set which is square summable. The deduction of this result assumes, as a sufficient (but not necessary) 
condition, that the mass function s{t) possesses a second derivative which is integrable in every compact subinterval 
of the time domain I. A time integration, over any such subinterval, of the partial sums of the square norms of the 
elements psp . combined with a suitable application of Luzin's theorem [31] (which is possible because the considered 
elements are measurable functions), shows then that the set formed by 2^ (namely, the ratios of the coefficients of /C.,-j) 
is square summable. Given the definition of Zf^ and relation pip , one can see that this implies that the set formed 
by the antilinear coefficients of K, A^;, is square summable as well 19]. But this last summability result is precisely 
the condition for the unitary implementability of the symplectomorphism K and, hence, of the unitary equivalence of 
the two Fock representations determined the complex structures J and Jq, related by that symplectomorphism. Since 
the discussion is valid for all admissible invariant complex structures J, we conclude that all such structures which 
besides allow for a unitary dynamics are indeed equivalent. Thus, our criterion of invariance under the three-torus 
isometrics and of unitarity in the evolution picks out a unique representation for the scalar field. 
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IV. UNIQUENESS OF THE FIELD DESCRIPTION 



We will show now that our criterion is not only capable of selecting a unique equivalence class of Fock representations 
for the quantization of a KG field with a time dependent mass in a flat spacetime with the spatial topology of a three- 
torus, but, beyond that, it also determines a unique canonical pair for the field among all those that are related 
by a linear canonical transformation varying in time and in which the field configuration gets scaled. This kind 
of scaling transformations are often found in cosmological contexts, either when dealing with test fields or with 
perturbations around homogeneous and isotropic solutions, which can represent a genuine background spacetime, an 
effective spacetime on which the propagation takes place once certain quantum effects are taken into account, or just 
an auxiliary background in which one formulates in a simpler form the field dynamics (for instance after dimensional 
reduction in symmetric models in General Relativity). The scaling absorbs part of the time dependence of the field, 
which is assigned to the time variation of the background. Although the classical formulations obtained for the field 
with this class of transformations are all equivalent, this ceases to be the case in the quantum theory, both because 
not all linear canonical transformations admit a unitary implementation and because the time dependence of the 
transformation changes the dynamics, in particular affecting its properties of unitarity. As a consequence, the criterion 
of a unitary evolution has different implications in the distinct formulations reached with these transformations. 



A. Time dependent canonical transformations 



The most general linear canonical transformation which includes a scaling of the field configuration and allows for 
a time dependence of the linear coefficients has the form 

= /(t)^, p^^^+g{t)^. (36) 

In the last equation, we have taken again into account that the determinant of the three-torus metric is equal to one. 
We assume that the real functions f{t) and g{t) that characterize the transformation are at least twice differentiable, 
in order to respect the differentiability properties in the field equations. Besides, we suppose that the function f{t) 
never vanishes, so that no spurious singularity is introduced in the field with the considered transformation. Finally, 
by means of a constant canonical linear transformation (which does not change the Fock representation of the field 
system), we can always set the initial values of the two time functions involved in our change equal to /(to) = 1 and 
9 (to) = [l6| . so that the original and the transformed canonical pairs coincide initially. 

For the new pair ((/>, P^), we start adopting the representation determined by the complex structure Jq. The time 
dependent transformation which leads to this pair changes the dynamics with respect to the original one, U. The new 



dynamical transformation U has 2x2 blocks labeled again by the LB eigenvalue number n, and given by [21[: 



Un{t,to)^C„{t)Un{t,to), (37) 

where 

r (f^ - f Mt) + ^Gnit) f-{t) + tGnit) \ 

2f±it):^fit)±^, G„(i):=f^. (39) 
fit) 2uj„ 

On the other hand, let us recall that the most general invariant complex structure J, compatible with the symplectic 
structure, has already been characterized in Sec. IIII Bl it is related with Jo by a symplectomorphism K of the type 
pop. In total, then, we see that the new dynamics U turns out to admit a unitary implementation with respect 
to an invariant complex structure J if and only if the beta Bogoliubov functions of the transformations with blocks 
lCZ^Cn{t)Un{t,to)ICfi are square summable over all the possible values of ft (the label of the sine and cosine modes). 
These beta functions, that we will call i3i{t,to), adopt an expression similar to that given in Eq. but with the 

Bogoliubov functions a„(t, to) and l3n{t, tg) -corresponding to the reference complex structure Jg- replaced with those 
of the evolution U for the pair {(f>,P^) (see Ref. (2l|): 

a„(t,to) - f+{t)an{t,to) + f-{t)p:,{t,to) + iGn{t)[an{t,ta) + I3:^{t,t„)], (40) 

Mt:to) = f+it)/3nit,to) + f-it)a*^{t,to)+iGnit)[[3nit,to)+alit,to)]. (41) 

In the rest of this section, we will demonstrate that a unitary dynamics is possible only if f{t) — 1 and g{t) = at all 
times, that is, if we describe our field precisely with the original canonical pair, associated with the formulation as a 
KG field in fiat spacetime with compact spatial sections and a time dependent mass. 
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B. Uniqueness of the scaling 



We will first show that f{t) must be the identity function in I. In order to do this, we simply adapt the proof 
explained in Ref. [2l| to the three-torus case. For each eigenvalue of the LB operator, —oj^, let us choose a value A/„ 
of the label ft among all those whose Euclidean norm as a vector coincides with a;„. We then consider the sequence 
with elements /^^- {t, to)- This sequence is a subset of the beta functions to), obtained by ignoring the degeneracy 
of the LB eigenspaces. Since the latter set is square summable if we admit that the dynamics is unitary, then, in 
particular, the sequence with labels M„ is also square summable at all times. Recalling the asymptotic behavior (fT2|) 
and employing that | k^^j | > 1 , it is not difficult to check that the considered square summability implies a vanishing 



limit at all times for the sequence with terms 



Ju>„{t-to) 



- 2iZj,j^ sin[w„(t - to)]f+{t). 



(42) 



We will now use that, for any given eigenvalue Wp, all its multiples by any positive integer I G N"*" again correspond 
to eigenvalues of the LB operator (since, given a mode with label n^, one can consider all modes with labels of the 
form Ifip). By restricting the possible values of a;„ in Eq. (|42l) to the subsequence lijjp, it is clear that one also gets a 



Let us call 



the value of z^g when w„ 



= lujp, and consider 



sequence with vanishing limit, this time when Z — > cxd. 

only subsets of instants of time of the form t = to + 2'Kr /ujp, with i e I and r any (allowed) fixed positive integer. 
From the real and imaginary parts of Eq. ()42[) , we then obtain that the sequences formed by the following terms must 
have a zero limit at infinitely large I: 



i-$n 



, , 2'Kr 
f-[to + 

LOp 



(43) 



The symbol denotes the real part. Notice that our conclusion about the limit of these terms is valid for every 
possible value of r and Wp. 

In Refs. , [m it is proven that, if the terms in Eq. (j42|) tend to zero, the two factors multiplying /_ [to + 2'Kr/ijjp) 
in the above equation cannot tend simultaneously to zero in any subsequence at large /. Therefore, we necessarily 
have that the value of /_ at to + ^-Krjujp is the number that must vanish. And this must be so for all possible values 
of the integer r and of the eigenvalue Wp. Since the set of points to + 2'Kr jujp obtained by varying those two quantities 
is dense in the considered interval I, and since /_ is continuous (because so is /, as well as 1//, for / cannot be zero), 
we conclude that it must coincide with the zero function in the whole of I. Then, recalling that jito) = 1, we get that 
/(f) has to be the unit function in its entire time domain. In this way, we reach the result that no scaling of the field 
configuration is allowed by our combined criterion of invariance under the three-torus isometries and the unitarity of 
the dynamics, as we wanted to prove. 



C. Uniqueness of the field momentum 



We will end this section by proving that no change in the momentum is permitted by the requirement of symmetry 
invariance and unitary evolution, so that the function g(t) must vanish. We return to the expression of the beta 
functions for the dynamics of the system after performing a time dependent linear canonical transformation, but now 
specialized to the case ][t) = 1, in accordance with the discussion of the previous subsection. The demand that, for 
an invariant complex structure J, the dynamics admit a unitary implementation amounts to the square summability 
of the set formed by iiiit^to) at all instants of time. Recalling that > 1, this summability ensures the same 
property for the set given by l3i{t,to) / {nVf' . Then, using the asymptotic relations (|12p . and that |z,-j| < 1, one can 
deduce that the set formed by 

Gn{t) {e»[""(*-*o)-*d + |^.|2g-.K(t-to)-5d + 2\zn\ cosK(t - to)]} + 2|z^jp K(t, to)] (44) 

is square summable at all times in the considered interval I. Here, 5fi is the phase of z^. Obviously, the square 
summability is also true for the set obtained by dividing those terms by a;„, since this eigenvalue provides a sequence 
that diverges to infinity. Hence, employing that gn/^^Jn is a summable sequence (as we proved in Sec. IIII Ap and the 
definition of Gn{t), we conclude that the set with elements jz^jp [a„(i, ^q)] /'^n niust also be square summable. By 
performing a convenient time average and making use of Luzin's theorem, one can show then that the set formed by 
\zn\/i-o„ has to be square summable as well. 
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Taking into account this result in the consideration of the terms (j44p . we arrive at the square summabihty of the 
set formed by 

G„(t)e^[""(*-*°)-^''l + 2\zn\3 K(t, to)] • (45) 
In particular, the imaginary part of these quantities, namely 

f^sin[cj„(t-to)-^d (46) 

[where we have used again the definition of G'„(t)], is also a square summable set at all times. If there existed a 
subinterval of I where the function g{t) did not vanish, a suitable time integration over it would lead to the conclusion 
that the sequence with elements gn/w^ must be summable [l^. But this sequence has in fact a divergent sum, because 
the sum exceeds that of Dn/{N + 1)^ over the positive integers, which clearly diverges, given the asymptotic behavior 
Dn oc N"^ discussed in Sec. IIII Al This eliminates the possibility that g{t) may differ from zero in any subinterval of 
L Since the function is continuous, this implies that g{t) has to be the zero function. 

Summarizing, our criterion determines completely the choice of canonical pair for the field among all the possibilities 
related by means of a time dependent linear canonical transformation. The criterion removes the freedom to scale 
the field, and to redefine its momentum by including a contribution that is linear in the field configuration. 



V. CONCLUSIONS 



We have analyzed two types of ambiguities that appear in the Fock quantization of scalar fields in cosmological 
spacetimes with spatial sections that are isomorphic to a three-torus. The first of these ambiguities is related to 
the possibility of scaling the field by means of a time dependent function, which assigns to the (physical, effective 
or auxiliary) background cosmological spacetime part of the time variation. This scaling can be viewed as resulting 
from a linear canonical transformation, in which the field momentum gets the inverse scaling (compared to the field 
configuration). In addition, the momentum may admit a linear contribution of the field configuration, with a time 
dependent coefficient. Canonical transformations of this kind often lead to a simpler and better behaved formulation 
for the system [5]. Each of these transformations provides a different canonical pair for the field description and, 
furthermore, changes the dynamics, since part of the time dependence is absorbed in the background. The other 
ambiguity that we have considered refers to the possible choices of Fock representation for the CCR's, once a specific 
canonical pair (and dynamical evolution) is given for the field. The physically different representations can be 
understood as corresponding to inequivalent choices of vacuum for the Fock construction. Alternatively, the physical 
freedom in the selection of a representation can be assigned to the possible inequivalent choices of complex structure. 
This ambiguity is well known in quantum field theory [2|, and it is common to remove it by introducing certain 
requirements on the vacuum, or equivalently on the complex structure, such as incorporating certain symmetries of 
the background spacetime, or presenting a especially good local or dynamical behavior. 

We have put forward a criterion to remove these ambiguities in systems that, by means of one of the considered 
time dependent canonical transformations, can be formulated as a KG field with time varying mass propagating in 
a fiat spacetime with (compact) three-torus spatial topology. The criterion consists of two requirements. First, that 
the complex structure (and hence the vacuum) be invariant under the isometrics of the three-torus, equipped with 
the standard metric. Second, that this complex structure also allows for the unitary implementation of the dynamics. 
Notice, in particular, that the considered ambiguity under time dependent canonical transformations affects the 
dynamical evolution of the system, so that our demand of unitarity has different implications for the distinct field 
descriptions obtained with such transformations. On the other hand, requiring a unitary dynamics guarantees that 
the standard probabilistic interpretation of quantum field theory is consistent in the evolution, so that one does not 
have to renounce to it. 

We have demonstrated that our criterion is indeed capable of removing the two mentioned types of ambiguities, 
selecting a unique canonical pair for the field among all those related by linear canonical transformations with time 
dependent coefficients and, furthermore, selecting a unique class of unitarily equivalent Fock representations for the 
corresponding CCR's. This class contains the representation that would be naturally associated with the case of a 
massless field, though now employed to define a Fock quantization of a system that is not only massive, but whose 
mass changes in time. To obtain this uniqueness result, we have assumed only a very mild restriction on the time 
dependent mass: that it must possess a second derivative which is integrable in every compact subinterval of the time 
domain. In fact, this assumption is just a sufficient condition, but does not even seem to be strictly a necessary one. 
Let us also emphasize that the conclusion about the uniqueness of the Fock quantization is valid for any possible time 
domain provided that it is a (non-infinitesimal) interval of the real line. 
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Our result about the uniqueness of the Fock quantization can be apphed to a number of physically interesting 
situations in cosmology. Current observations indicate that the large scale structure of the universe is (approximately) 
homogeneous and isotropic. Therefore, (quantum) matter fields in cosmology are naturally described by quantum 
field theory in FRW spacetimes. Besides, inhomogeneities can then be treated as perturbations in a fairly good 
approximation At leading order, these perturbations are also described as linear fields that propagate in FRW 
spacetimes. Furthermore, the observations favor a spatially flat cosmology, which is precisely the case studied in 
this work. On the other hand, the fact that the spatial topology is assumed to be compact should not pose a 
severe restriction, because one would expect that, beyond a certain cosmological scale related with the Hubble radius, 
physical interactions should have no relevant effect. Then, the physics would not be altered importantly by considering 
a compactification scale larger than this cosmological one. In this general context, a simple system that can be 
reformulated by means of a scaling as a KG field with time dependent mass in a flat spacetime is a massive, minimally 
coupled scalar field. A description of this type is also obtained after a suitable scaling for cosmological perturbations in 
conformal time. For instance, this is the case of tensor perturbations (for which one can generalize our considerations, 
presented here for scalar fields) and of the gauge invariant energy density perturbation amplitude 0, [s^] ■ In addition, 
one obtains a similar kind of field description for scalar perturbations of a massive scalar field around flat FRW 
spacetimes after adopting (e.g.) a longitudinal gauge, with the only caveat that the KG equation is modified with 
subdominant terms which, nonetheless, do not affect the asymptotic behavior employed in our discussion @, H^l- The 
fact that, in all these cases, the criterion of spatial symmetry invariance and unitary evolution turns out to determine 
a unique Fock quantum description provides the quantum field theory and its predictions with the desired robustness. 
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